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Abstract. We present new results on the cosmological evolution of unabsorbed (type-1) active galactic nuclei 
(AGN) selected in the soft (0.5-2 keV) X-ray band. From a variety of ROSAT, XMM-Newton and Chandra 
surveys we selected a total of ~ 1000 AGN with an unprecedented spectroscopic and photometric optical/NIR 
identification completeness. For the first time we are able to derive reliable space densities for low-luminosity 
(Seyfert-type) X-ray sources at cosmological redshifts. The evolutionary behaviour of AGN shows a strong de- 
pendence on X-ray luminosity: while the space density of high-luminosity AGN reaches a peak around z ~ 2, 
similar to that of optically selected QSO, the space density of low-luminosity AGNs peaks at redshifts below 
z — 1. This confirms previous ROSAT findings of a luminosity-dependent density evolution. Using a rigorous 
treatment of the optical identification completeness we are able to show that the space density of AGN with 
X-ray luminosities L x < 10 45 erg s _1 declines significantly towards high redshifts. 
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1. Introduction 

In recent years the bulk of the extragalactic X-ray 
background in the 0.1-10 keV band has been resolved 
into discrete sources with the deepest ROSAT, Chandra 
and XMM-Newton observations (Hasinger et al. [1998 
Mushotzky et al. 120001 Giacconi et al. l200Tl 120021 Hasinger 
et al. I200T1 Alexander et al. 120051 Worsley et al. 12004) . 
Optical identification programmes with Keck (Schmidt et 
al. HOOSl Lehmann et al. I200T1 Barger et al. I200T1 I2003}l 
and VLT (Szokoly et al. 12001 Fiore et al. l2003jl find pre- 
dominantly unobscured AGN-1 at X-ray fluxes Sx > 
and a mixture of unobscured AGN- 



10 14 erg cm 2 s 

1 and obscured AGN-2 at fluxes 10" 14 > S x > 10 



-15.5 



erg cm s with ever fainter and redder optical coun- 
terparts, while at even lower X-ray fluxes a new popula- 
tion of star forming galaxies emerges (Hornschcmcicr ct al. 
120051 Rosati et al. 120021 Norman et al. 15003} . At optical 
magnitudes R>24 these surveys suffer from large spectro- 
scopic incompleteness, but deep optical/NIR photometry 
can improve the identification completeness significantly, 
even for the faintest optical counterparts (Zheng et al. 
Mainieri et al. 
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The AGN/QSO luminosity function and its evolution 
with cosmic time are key observational quantities for un- 
derstanding the origin of and accretion history onto su- 
permassive black holes, which are now believed to occupy 
the centers of most galaxies. X-ray surveys are practi- 
cally the most efficient means of finding active galactic nu- 
clei (AGNs) over a wide range of luminosity and redshift. 
Enormous efforts have been made by several groups to fol- 
low up X-ray sources with major optical telescopes around 
the globe, so that now we have fairly complete samples of 
X-ray selected AGNs. In this work we concentrate on un- 
absorbed (type-1) AGN selected in the soft (0.5-2 keV) 
X-ray band, where due to the previous ROSAT work (see 
Miyaji et al.. 120001 120PT1 hereafter Paper I and II) com- 
plete samples exist, with sensitivity limits varying over 
five orders of magnitude in flux, and survey solid angles 
ranging from the whole high galactic latitude sky to the 
deepest pencil-beam fields. These samples enable us to 
construct and probe luminosity functions over cosmolog- 
ical timescales, with an unprecedented accuracy and pa- 
rameter space. 

Conceptually, space densities and luminosity functions 
are simply derived by dividing the observed number of ob- 
jects N by the volume V, in which they have been sur- 
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veyed. The binned luminosity function derived in this fash- 
ion generally does not represent the center of the bin. In 
paper II, we introduced the estimator N ^, s /N mt ^, where 
•^mdi is the number of objects expected from an analyti- 
cal representation of the luminosity function. Scaling the 
model value of the analytical function at the bin center 
by the estimator removes the binning bias to first order. 
This method is applied in Sect. 0] A quite different method 
based on 1/V max values (Schmidt .1968;) for individual ob- 
jects is used in Sect. |S] It involves a derivation of the 
zero-redshift luminosity function that is free of binning 
bias. The luminosity function at higher redshifts is again 
derived by employing an analytical representation of the 
luminosity function and scaling it by the ratio of observed 
over expected numbers. The use of individual 1/V max val- 
ues in this case allows accounting for an effective optical 
magnitude limit beyond which redshifts have generally not 
been obtained in some surveys. 

Throughout this work we use a Hubble constant 
Hq = 70/i7o kms^Mpc -1 and cosmological parameters 
(Q m , = (0.3, 0.7) consistent with the WMAP cosmol- 
ogy (Spergel et al-EESJ- 

2. The X-ray selected AGN-1 sample 

For the derivation of the X-ray luminosity function and 
cosmological evolution of AGN we have chosen well- 
defined flux-limited samples of active galactic nuclei, with 
flux limits and survey solid angles ranging over five and six 
orders of magnitude, respectively. To be able to utilize the 
massive amount of optical identification work performed 
previously on a large number of shallow to deep ROSAT 
surveys, we restricted the analysis to samples selected in 
the 0.5-2 keV band. In addition to the ROSAT surveys 
already used in Paper I and II, we included data from 
the recently published ROSAT North Ecliptic Pole Survey 
(NEPS, Gioia et al. 120351 Mullis et al. IMHjl . from an 
XMM-Newton observation of the Lockman Hole (Mainieri 
et al. 1271172) and the Chandra Deep Fields South (CDF-S, 
Szokoly et al. 127701 Zheng et al. 127701 Mainieri et al. lMTS) 
and North (CDF-N, Barger et al. l2I70"Tll2UU51l . In order to 
avoid systematic uncertainties introduced by the varying 
and a priori unknown AGN absorption column densities 
we selected only unabsorbed (type-1) AGN, classified by 
optical and/or X-ray methods. We are using here a defini- 
tion of type-1 AGN, which is largely based on the presence 
of broad Balmer emission lines and small Balmer decre- 
ment in the optical spectrum of the source (optical type-1 
AGN, e.g. the ID classes a, b, and partly c in Schmidt 
et al. 1199811 . which largely overlaps the class of X-ray 
type-1 AGN defined by their X-ray luminosity and un- 
absorbed X-ray spectrum (Szokoly et al. I2UU4|I . However, 
as Szokoly et al show, at low X-ray luminosities and in- 
termediate redshifts the optical AGN classification often 
breaks down because of the dilution of the AGN excess 
light by the stars in the host galaxy (see e.g. Moran et 
al. lj2D02|l . so that only an X-ray classification scheme can 
be utilized. Schmidt et al. ( 1993J) have already introduced 



the X-ray luminosity in their classification. For the deep 
XMM-Newton and Chandra surveys we in addition use 
the X-ray hardness ratio to discriminate between X-ray 
type-1 and type-2 AGN, following Szokoly et al. (2004). 

In order to convert the count rates observed in the 0.5- 
2 keV band to unabsorbed 0.5-2 keV fluxes, we assumed a 
power law AGN spectrum with a photon index of T = 2.0 
and Galactic absorption. Typical values of the Galactic 
neutral hydrogen column density are (0.5 — 1) x 10 20 cm~ 2 
for the deep surveys and a maximum of 16 x 10 20 cm -2 for 
a small portion of the sky covered by the ROSAT Bright 
Survey (RBS). Because the band, in which the AGN are 
selected, is the same as the one for which we calculate 
the fluxes, systematic differences in the true AGN power 
law indices have a negligible effect on the derived fluxes. 
Assuming spectral indices in the range T = 1 — 3, the 
conversion between the observed 0.5-2 keV count rates 
and the X-ray flux S x (here and hereafter, S x represents 
the 0.5-2 keV flux and 5*xi4 is the same quantity measured 
in units of 10~ ergs -1 cm -2 ) varies by less than 10%. 

The surveys we have used are summarized in Tabled 
A total of 944 X-ray selected type-1 AGN were compiled 
from eight independent samples containing a total of 2566 
soft X-ray sources. The number of unidentified sources 1 
in these samples is only 86 (of which 57 could be AGN-1), 
yielding an unprecedented identification fraction of 97%. 
Due to the extreme faintness of the optical counterparts, 
the lowest identification fractions are achieved in the re- 
cent deepest samples: 87% for the XMM-Newton survey 
in the Lockman Hole and 88% in the CDF-N. A surpris- 
ingly high identification fraction of 98% has been achieved 
in the CDF-S through the utilization of photometric red- 
shifts based on extremely faint optical/NIR photometry. 

For the computation of the soft X-ray luminosity func- 
tion SXLF, it is important to define the available survey 
solid angle as a function of limiting flux. In case there is 
incompleteness in the spectroscopic identifications in the 
ROSAT surveys, we have made the usual assumption that 
the redshift/classification distribution of these unidenti- 
fied sources is the same as the identified sources at simi- 
lar fluxes by defining the 'effective' survey solid angle as 
the geometrical survey solid angle multiplied by the com- 
pleteness of the identifications (see Paper I). This assump- 
tion is not correct when the source is unidentified due to 
non-random causes, in particular its optical faintness. The 
treatment of this identification incompleteness and the ef- 
fect of the optical limit on the derived space densities is 
discussed in detail in Sect El 

In several surveys we had to choose an X-ray flux 
limit a posteriori, based on optical completeness criteria, 
i.e. maximising the number of optically identified sources, 
while simultaneously minimising the number of unidenti- 
fied objects. This procedure can introduce a bias against 

1 We call unidentified sources those, which do not have a 
reliable redshift determination, either through spectroscopy or 
through photometric redshifts; however practically all of these 
have optical or NIR counterparts 
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Table 1. The soft X-ray sample 



Survey a 


Solid Angle 


S'xi4,lim 


iVtot 


A'agn— i b 


M ' umd 




[deg 2 ] 


[cgs] 








RBS 


20391 


» 250 


901 


203 





SA-N 


684.0-36.0 


47.4-13.0 


380 


134 


5 


NEPS 


80.8-70.5 


12.4-10.1 


252 


101 


1 


RIXOS 


19.5-15.0 


10.2-3.0 


340 


194 


14 


RMS 


0.74-0.32 


1.0-0.5 


124 


84 


7 


RDS/XMM 


0.126-0.087 


0.38-0.13 


81 


48 


8 


CDF-S 


0.087-0.023 


0.022-0.0053 


293 


113 


1 


CDF-N 


0.048-0.0064 


0.030-0.0046 


195 


67 


21 


Total 






2566 


944 


57 



a Abbreviations - RBS: The ROSAT Bright Survey (Schwope et al.|2000); SA-N: ROSAT Selected Areas North (Appenzeller et 
al. lT9T38l : NEPS: ROSAT North Ecliptic Pole Survey (Gioia et al.. l3D03t : RIXOS: ROSAT International X-ray Optical Survey 
(Mason et al. l2000ll . RMS: R OSAT Medium Deep Survey, consisting of deep PSPC pointings at the North Ecliptic Pole (Bower 
et al. 19961, the UK Deep Survey (McHardy et al. 11998^ . the Marano field (Zamorani et al. I1999(l and the outer parts of the 
Lockman Hole (Schmidt et al. 119981 Lehmann et al. 1200011 : RDS/XMM: ROSAT Deep Survey in the central part of the Lockman 
Hole, observed with XMM-Newton (Lehmann et al. 120011 Mainieri et al. 120021 Fadda et al. 1200211 : CDF-S: The Chandra Deep 
Field South (Szokoly et aL EUOH Zheng et al. 1277011 Mainieri et al. I30"0"5l : CDF-N: The Chandra Deep Field North (Barger et 

ai-nssunssni- 

b Excluding AGNs with z < 0.015. 

c Objects without redshifts, but hardness ratios consistent with type-1 AGN. 



optically faint sources, if the reason for the missing red- 
shift is the optical faintness of the source and in fields 
with relatively few objects. We have tried to minimise the 
impact of this "gerrymandering" effect, e.g. by allowing 
a number of unidentified sources to enter the sample and 
then defining the corresponding X-ray flux limits in the 
geometric mean between the last identified and the next 
unidentified source. In addition, the wide range for X-ray 
and thus optical flux limits in our survey tends to reduce 
biases, which occur at the flux limit of individual surveys. 
The problem of missing redshifts in the faintest surveys is 
addressed specifically in Sect 0] and [5] 

Below we summarize our sample selection and com- 
pleteness for each survey. Figurenshows the AGN-1 sam- 
ple in the redshift - luminosity plane. Figure |3 gives the 
combined solid angle versus flux curve. Both the sample 
and the solid angle coverage are available in computer 
readable form under http : //mpe .mpg . de/~ ghasinger. 

2.1. The ROSAT Bright Survey (RBS) 

The RBS identified the brightest ~ 2000 X-ray sources 
detected in the ROSAT All-Sky Survey (RASS, Voges 
et al. I1999[l at high galactic latitude, |6| > 30°, exclud- 
ing the Magellanic Clouds and the Virgo cluster, with 
ROSAT PSPC count rates above 0.2s" 1 . This program 
achieved a spectroscopic completeness of 99.5% (Schwope 
et al. I2000[) . We selected the sub-sample of 931 sources 
with count rates above 0.2 s" 1 in the ROSAT 0.5-2 keV 
band (PSPC channels 52-201), which is 100% identified. 
Since the absorption in our galaxy varies from place to 
place, the same count rate limit corresponds to different 
0.5-2 keV flux limits based on the different galactic Nn 
values. The Nn value ranges from (0.5 — 16) x 10 20 cm~ 2 
in the RBS survey area. Correspondingly, the survey solid 



angle varies steeply with flux from about 3000 deg 2 at a 
flux limit of 5x14 = 246 to a total of 20391 deg 2 at a flux 
limit of Sxu — 360. 

2.2. The RASS Selected-Area Survey North (SA-N) 

This survey gives optical identifications of a represen- 
tative sample of northern (<5 > —9°) RASS sources in 
six study areas outside the Galactic plane (\b\ > 19.6°) 
with a total of 685 deg 2 . A count rate limited complete 
RASS subsample comprising 674 sources has been identi- 
fied (Appenzeller et al. I1998f) . The fields selected for the 
survey have a Galactic column density in the range Nn = 
(2— 11) x 10 20 cm" 2 . We have further selected our sample 
such that each of the six fields has a complete ROSAT 
hard-band (0.5-2 keV; channels 52-201) countrate-limited 
sample with complete identifications (Ci?o.5-2koV > 0.01- 
0.05 cts s" 1 ). To avoid overlap with the RBS, those sources 
common in both samples were removed from SA-N, yield- 
ing a total of 406 sources with 98.5% spectroscopic com- 
pleteness. 

2.3. The ROSAT North Ecliptic Pole Survey (NEPS) 

The RASS data in a contiguous area of 80.7 deg 2 around 
the North Ecliptic Pole (Galactic latitude b > 29.8°) have 
been used to construct a survey consisting of 445 X-ray 
sources detected above a 4cr threshold. Gioia et al. ( 2003) 
and Mullis et al. ijMlljl have identified 99.6% of these 
sources and determined redshifts for the extragalactic ob- 
jects. Since the exposure in the ROSAT All-Sky Survey in- 
creases significantly towards the North-Ecliptic Pole, the 
actual survey sensitivity is a strong function of ecliptic 
latitude. The original NEPS sample is selected in the full 
ROSAT PSPC band. For consistency with the other sur- 
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Fig. 1. The AGN-1 soft X-ray sample in the 2-log L x plane. 
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Fig. 2. The survey solid angle of the combined soft X-ray sample as a function of flux. The solid line shows the case 
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veys used in our work, we selected sources detected sig- 
nificantly in the PSPC hard band (0.5-2 keV; channels 
52-201) by specifying a hard count rate limit as a func- 
tion of ecliptic latitude. New 0.5-2 keV fluxes were calcu- 
lated from the hard PSPC count rates in the same way 
as for the RBS and SA-N objects, taking into account 
the Galactic neutral hydrogen column density varying in 
the range 2.6 — 6.2 x 10 20 cm~ 2 across the NEPS region. 
Due to the large gradient in exposure times, we have cut 
the sample to 252 sources with fluxes above 5xi4 = 10.1, 
where the solid angle of this survey is 70.5 deg 2 , increas- 
ing to 80.8 deg 2 at 5xi4 = 12.4. Only one of these sources 
remains unidentified in the NEPS sample. 

2.4. The ROSAT International X-ray/Optical Survey 
(RIXOS) 

The ROSAT International X-ray/Optical Survey 
(RIXOS, Mason et al. I2000J) is a medium-sensitivity 
survey and optical identification program of X-ray 
sources discovered in R OS A Thigh Galactic latitude fields 
(|6| > 28 deg) and observed with the Position Sensitive 
Proportional Counter (PSPC) detector with a minimum 
exposure time of 8 ks. The survey comprises 82 ROSAT 
PSPC fields and made use of the central 17 arcmin 
of each field, however, excluding the target region for 
pointings on known X-ray sources. The total survey 
contains 395 X-ray sources, selected in the PSPC 0.5-2 
keV band. A flux limit of 5xi4 = 3.0 was adopted for 
the survey, substantially above the detection threshold of 
each field, however, the actual spectroscopic completeness 
limit varies from field to field. We have chosen a strategy, 
which on one hand maximises the sample of identified 
AGN-1 and on the other hand minimises the number 
of unidentified sources. There are 51 fields (12.3 deg 2 ) 
identified completely down to the survey flux limit. Three 
fields have such a low identification fraction, that we 
ignore them. For the remaining 28 fields we allow at most 
one unidentified source. If an unidentified source has the 
lowest flux of the subsample of a particular field, we 
exclude this source and raise the flux limit for this field to 
the geometric average between the flux of this source and 
that of the last identified source. This way we can define 
a clean RIXOS sample comprising 340 objects and only 
14 unidentified sources, i.e. an identification fraction of 
95.9%. The survey solid angle, corrected for spectroscopic 
incompleteness, rises from 15.0 deg 2 at 5xi4 = 3.0 to 19.5 
deg 2 at 5xi4 = 10.2. 

2.5. The ROSAT Medium Survey (RMS) 

For this work we have grouped a number of medium-deep 
i? 05^4 X" surveys with flux limits in the range 5x14 = 0.5—1 
into the RMS. In particular these comprise pointed obser- 
vations at the North Ecliptic Pole (Bower et al. Q996), 
the UK Deep Survey (McHardy et al. lT^5|l . the Marano 
field (Zamorani et al. I1999f) and the outer parts of the 



Lockman Hole (Schmidt et al. 119981 Lehmann et al.EDOO). 
The North Ecliptic pole pointing covers the same sky area 
as the center of the NEPS, however, to a flux limit of 
S"xi4 = 1-0. Again, we remove the overlapping sources be- 
tween the two surveys. For the UK Deep Survey and the 
Marano Field we define a flux limit of 5xi4 = 0.5, follow- 
ing Paper I. For the ROSAT PSPC survey of the Lockman 
Hole we only chose the region not covered by the deeper 
RDS/XMM survey (see section^))) but otherwise selected 
the completely identified sample with the same flux lim- 
its as those chosen for the ROSAT Ultradeep Survey UDS 
by Lehmann et al. IWiTI 5x14 > 0.96 for PSPC off-axis 
angles in the range 12.5-18.5 arcmin and 5xi4 > 0.55 
for off-axis angles smaller than 12.5 arcmin. Overall, the 
RMS contains 124 sources, at an identification complete- 
ness of 94.4%. Correspondingly, the corrected survey solid 
angle varies in the range 0.30-0.70 deg 2 for flux limits 
5xi4=0.5-l. 

2.6. Deep XMM-Newton survey of the Lockman Hole 
(XMM/RDS) 

The Lockmam Hole (XMM/RDS) has been observed by 
XMM-Newton a total of 17 times during the PV, AO-1 
and AO-2 phases of the mission, with total good expo- 
sure times in the range 680-880 ks in the PN and MOS 
instruments (see Hasinger et al. l2001ll2004l and Worsley et 
al. 120041 for details). Spectroscopic optical identifications 
of the R OS A T sources in the LH have been presented by 
Schmidt et al. (fTM^ and Lehmann et al. lj2U0Tjl EMTjl 
and a new catalogue from the XMM-Newton PV phase 
is given in Mainieri et al. (|2002p . Some photometric red- 
shifts have been discussed in Fadda et al., H2002|l . Here 
we selected sources from the 770 ksec dataset (Brunner et 
al., 2005, in prep.) with additional spectroscopic identifi- 
cations obtained with the DEIMOS spectrograph on the 
Keck telescope in spring 2003 and 2004 by M. Schmidt 
and P. Henry (Szokoly al., 2005, in prep.). In order to 
maximise the spectroscopic/photometric completeness of 
the sample, we selected objects in two off-axis intervals: 
5xi4 = 0.38 for off-axis angles in the range 10.0-12.5 ar- 
cmin and 5xi4 = 0.13 for off-axis angles smaller than 10.0 
arcmin. The total number of sources in the XMM/RDS is 
81, with 8 potential AGN-1 still unidentified. 

2.7. Chandra Deep Field South (CDF-S) 

We have used the catalogue of Giacconi et al. (2002) based 
on the 1 Ms observation of the CDF-S (Rosati et al. 
2002]). Spectroscopic identifications with the FORS instru- 
ments at the ESO VLT have been obtained by Szokoly et 
al. (2004), yielding a spectroscopic completeness around 
60%. Additional spectroscopic redshifts of CDF-S X-ray 
sources have been obtained with the VIMOS spectrograph 
at the ESO VLT (Lefevre et al., 2004). The field is also in- 
cluded in the COMBO- 17 intermediate-band optical sur- 
vey, which gives very reliable photometric redshifts for the 
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brighter sources (Wolf et al., I2004fl . Very deep NIR pho- 
tometry has been obtained with the ISAAC camera at 
the VLT in conjunction with deep optical imaging with 
the HST ACS as part of the GOODS project (Dickinson 
& Giavalisco E001 Mobasher et al. rMHf) . The CDF-S 
therefore offers the highest quality photometric redshifts 
of faint X-ray sources, which are discussed in Zheng et 
al. ifSfMj) and Mainieri et al. p¥)5)l . Using the ISAAC 
images, tentative photometric redshifts could even be as- 
signed to several of the extreme X-ray /optical sources 
(EXOs) discussed by Koekemoer et al. H2004[) . We selected 
all sources from the Giacconi et al. 120021 catalogue within 
10 arcmin from the Chandra pointing center significantly 
detected in the 0.5-2 keV band. The sample thus contains 
a total of 293 objects. Combining all spectroscopic and 
photometric redshifts, only 2 sources in the CDFS remain 
unidentified, of which one could be an AGN-1. The sur- 
vey solid angle for the CDF-S has been estimated using 
a simple off-axis dependent flux limit. The solid angle in- 
creases from 0.023 deg 2 at S Xl4 = 0.0053 to 0.087 deg 2 at 
Sxi4 = 0.027. 

2.8. Chandra Deep Field North (CDF-N) 

We have used selected X-ray sources from the 2 Ms CDF- 
N source catalogue by Alexander et al. (2003| along with 
optical identifications by Barger et al. J2003) for our 
AGN-1 sample. Following Szokoly et al. I|2004fl . we se- 
lected AGN-1 either from broad permitted Balmer lines or 
from the X-ray luminosity and hardness, using HR<-0.2. 
We set our flux limits such that we also have sufficient hard 
(2-8 keV) sensitivity to exclude objects with HR>-0.2 and 
to include as many sources as possible which meet these 
criteria. Unlike the CDF-S, the CDF-N exposure map has 
a complicated structure and a simple off-axis dependence 
of the limiting flux is not a good approximation. Thus we 
have used the rectangular region of 170 arcmin 2 , which 
has the deepest coverage and is mostly co-spatial with 
the region covered by HST ACS in the GOODS project 
(Giavalisco et al. I2UU1 Cowie et al. I2U03) |. Within this 
region, the exposure and background are smooth enough 
that the photon counts limit of the detected sources can 
be approximated by a simple function of off-axis angle. 
In practice, due to statistical fluctuations, three sources 
have upper limits to HR between -0.1 and -0.2 and those 
have been considered to meet our hardness ratio crite- 
rion. Among the 128 sources meeting the soft counts limit 
and hardness ratio criteria, 20 are unidentified and 5 are 
stars (85% completeness). Only one broad-line AGN had 
a harder hardness ratio than our limit; this was also in- 
cluded in our type-1 AGN sample. The flux-solid angle 
relation has been calculated from the "limiting flux map" , 
where the counts limit is divided by the soft-band ex- 
posure map (in seconds) and multiplied by the conver- 
sion factor of 5 • 10 -12 ergs -1 cm -2 (Alexander et al., 
2003 ) . Due to the incompleteness in this field, where most 
unidentified sources are optically faint, this sample has 



not been included in the analysis in Section^] but consid- 
ered in Sectional where a method is developed to account 
for the optical magnitude limit in calculating the survey 
volume. 

3. Number counts for different source classes 

The combination of a large number of surveys with a wide 
range of sensitivity limits and solid angle coverage presents 
a unique resource. On one hand, the surveys presented 
here resolve the soft X-ray background almost completely. 
On the other hand, we have an almost complete opti- 
cal identification and redshift determination for all con- 
stituents. For the first time in any astrophysical waveband 
we are thus in the position to study the complete contri- 
bution of different object classes to the X-ray background 
and their evolution with cosmic time. Using the solid an- 
gle versus flux limit curve given in Figure [5] we compiled 
number counts for the total sample including all classes 
of sources and for the subclass of AGN-1. Figure (top) 
shows the cumulative source counts. For clarity we also 
show normalized differential source counts dA^ci<Sxi4>Sxi4 
in the bottom panel of Figure 01 Euclidean source counts 
would correspond to horizontal lines in this graph. 

For the total source counts, the well-known broken 
power law behaviour is confirmed with high precision. We 
fitted a broken power law to the differential source counts 
and obtain power law indices of aj, — 2.34 ± 0.01 and 
a.f = 1.55±0.04 for the bright and faint end, respectively, 
a break flux of Sxi4 = 0.65 ± 0.10 and a normalisation 
of dN/dSxi4 = 103.5 ± 5.3 deg" 2 at S X u = 1-0 with 
a reduced x 2 = 1.51. The total differential source counts, 
normalized to a Euclidean behaviour (dN/dSxi4 x S* x f 4 
is shown with open symbols in Figure E3 We see that the 
total source counts at bright fluxes, as determined by the 
ROSAT All-Sky Survey data, are significantly flatter than 
Euclidean, consistent with the discussion in Hasinger et al. 
IPM^I . Moretti et al. (|2003j) . on the other hand, have de- 
rived a significantly steeper bright flux slope (a.b ~ 2.8) 
from ROSAT HRI pointed observations. This discrepancy 
can probably be attributed to the selection bias against 
bright sources, when using pointed observations where the 
target area has to be excised. 

Type-1 AGN are the most abundant population of soft 
X-ray sources. For the determination of the AGN-1 num- 
ber counts we include those unidentified sources, which 
have hardness ratios consistent with AGN-1 (a contribu- 
tion of ~ 6%, see Tableful. Figure shows, that the break 
in the total source counts at intermediate fluxes is pro- 
duced by type-1 AGN, which are the dominant popula- 
tion there. Both at bright fluxes and at the faintest fluxes, 
type-1 AGN contribute about 30% of the X-ray source 
population. At bright fluxes, they have to share with clus- 
ters, stars and BL-Lac objects, at faint fluxes they com- 
pete with type-2 AGN and normal galaxies. We fitted a 
broken power law to the differential AGN-1 source counts 
and obtain power law indices of a\, = 2.55 ± 0.02 and 
aj = 1.15±0.05 for the bright and faint end, respectively, 
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Fig. 3. (a) Cumulative number counts N(>S) for the total sample (upper thin line) and the AGN-1 subsample 
(lower thick line), (b) Differential number counts of the total sample of X-ray sources (open squares) and the AGN-1 
subsample (filled squares). The dot-dashed lines refer to broken powerlaw fits to the differential source counts (see 
text). The dashed red line shows the prediction for type-1 AGN based on the model described in in sectional 



a break flux of 5xi4 = 0.53 ±0.05, consistent with that of 4. The SXLF and the Space Density Function 

the total source counts within errors, and a normalisation 
of of dN/dSxu = 83.2 ± 5.5 deg~ 2 at S X i4 = 1.0 with a 
reduced x 2 = l-26. The AGN-1 differential source counts, 
normalized to a Euclidean behaviour (dN/dSxi4 x S^'ik) 
is shown with filled symbols in Figure [21 Also shown are 
the predictions of the best-fit SXLF models discussed in 
Section |S1 



4.1. Basic method 

In this section, we present the binned Soft X-ray 
Luminosity Function (SXLF) of type-1 AGNs. The ba- 
sic approach is to use the N° hs /N mdl estimator described 
in Paper II. The procedure is outlined below: 



1. Divide the combined sample into several redshift shells. 
For each redshift shell, fit the AGN XLF with a smooth 
analytical function using a Maximum-likelihood fit 
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over each object (i.e., without binning; see Paper I for 
details). 

2. For the fitted model in each redshift shell, check the ab- 
solute goodness of fit with one- and two-dimensional 
Kolmogorov-Smirnov tests (hereafter, 1D-KS and 2D- 
KS tests respectively; Press et al. 119921 Fassano & 
Franceschini I1987j) . The K-S tests are also for the un- 
binned data sets and thus are free from artefacts and 
biases due to binning. 

3. For each redshift shell, bin the objects in luminos- 
ity bins to determine the observed number of objects 
(N ohs ). 

4. For each luminosity bin, evaluate the analytical fit at 
the central luminosity /redshift (d$ mdl /cflogL x ). 

5. Calculate the predicted number of AGNs in the bin 
(N mdl ). 

6. The final result is 



d$/dlogi x = d$ mdl /d log£ x • N ODS /N 



obs , 



r mdl 



(i) 



For the analytical expression of the SXLF in each red- 
shift shell, we use the smoothed two power law formula. 
Because the redshift shells have a finite widths, the fit 
results depend on the evolution of the SXLF within them: 



d*(L x ,z) 
d log L x 



72-1 



e d (z,L x ), (2) 



where ed(z,L x ) is the density evolution factor. While 
the final results are insensitive to the detailed behav- 
ior of e&{z, L x ) within the shell at most locations in the 
(L x , z) space, we have taken our best-estimate by us- 
ing the luminosity-dependent density evolution (LDDE) 
model derived later in Sect. 14.41 The luminosity range of 
the fit is from logL x =42.0 to the maximum available lu- 
minosity in the sample. 

In this section, we tried to make the sample as com- 
plete as possible, and we excluded the CDF-N from the 
analysis, where the incompleteness fraction is significant 
and most of the unidentified sources are optically faint. 
All of the unidentified sources in the ROSAT samples are 
optically bright and the reasons for them to be unidenti- 
fied are mostly by random causes, i.e., are not correlated 
with the intrinsic properties of the source. For the CDF-S, 
extensive photometric redshift studies including COMBO- 
17 (Wolf et al.. l2004f> and a careful individual photometric 
redshift determination of X-ray sources by Zheng et al. 
2004 and Mainieri et al.. 120051 has left only one potential 
AGN-1 without redshift information. For RDS/XMM, 2 
of the 8 unidentified sources which could be type-1 AGNs 
from X-ray hardness/spectra criteria are optically bright 
and they have remained unidentified so far for random rea- 
sons. The remaining 6 are optically faint (i? > 24.0) and 
the reason for remaining unidentified may well be corre- 
lated with redshift. 

As our nominal case, we took the first-order approach 
and defined "effective" survey solid angle (as a function 
of flux) , which is the geometrical survey solid angle mul- 
tiplied by the completeness, i.e. the fraction of identified 



X-ray sources in the survey, whether or not they are op- 
tically faint or bright. The correction has been made in 
each survey. In addition, we have also considered the up- 
per bounds on the binned SXLF and the space density 
function from the sample where all the unidentified op- 
tically faint (i? > 24.0) sources in turn are assigned the 
central redshift of each redshift shell. 

In the latter case, we used the geometrical solid angle 
for the CDF-S and the incompleteness correction to the 
RDS/XMM solid angle was only for the optically bright 
R < 24 unidentified sources. 

4.2. The binned SXLF 

The best fit parameters of Eq. [21 for each redshift shell are 
shown in Table[3along with the results of the ID- and 2D- 
KS tests (see the notes of the table). The normalization is 
defined by: 



d $ (L x = 10 44 erg s 



Zc) 



d log L x 



(3) 



where z c is the central redshift of the shell. The parame- 
ter errors in Table [5] correspond to a likelihood change of 
2.7 (90% confidence errors), except for the normalization 
A44, which cannot be a fit parameter in the maximum- 
likelihood method. The errors of A 44 are simply taken 
as the 90% Poisson errors of the number of the sources. 
Defining the normalization at a fixed luminosity (log L K = 
44) minimizes its dependence on other parameters. In any 
case, Eq. |2 gives a statistically satisfactory expression for 
all redshift shells as shown in Table 

We have made luminosity bins starting with a min- 
imum luminosity of log L x = 42.0 with a smallest bin 
size of A log L x = 0.25 in each redshift shell. If there are 
fewer than 10 AGNs in a bin, we have further rebinned up 
to a maximum bin size of Alog L x — 1.0. Table [3] shows 
the full binned results for the nominal case, along with ob- 
served number of AGNs (A^ obs ), model (Tabled predicted 
number (N 111 ^) and final estimated values SXLF value at 
the center of each bin. For reference, the additional num- 
ber of AGNs for the case where all the optically-faint 
unidentified sources are assigned the central redshift of the 
bin (in duplicate, as described above) are also shown in 
the last column of Table The full SXLF in the 6 red- 
shift shells is plotted in Figs. in separate panels. In all 
but the closest redshift shell panel, the best-fit two power 
law function to the 0.015 < z < 0.2 SXLF (Table 0) are 
also overplotted for reference. Three overall analytical ex- 
pressions discussed in Sect. 14.41 are also overplotted for 
comparison as discussed there. Because of the high com- 
pleteness of our sample, the redshift distribution of the 
optically faint sources affects the final SXLF results very 
little except for the case where all of them happen to fall 
into the highest redshift shell. In this case, the SXLF in 
the 3.2 < z < 4.8, 44 < log L x bin almost double. This 
is also verified by a comparison with the alternative ap- 
proach outlined in Sect.[5l 
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Table 2. Best-fit parameters for each redshift shell 
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z-range z~ c N A% 4 log LJ,* 7i 72 KS-prob b 



0.0-0.2 


0.1 
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■ io- 


-7 


43.45+af 7 
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o 1+0.4 
z - x -0. 3 


0.84,0.17,0.14 


0.2-0.4 


0.3 


139 


(9.76 ±0.83) 


10" 


-7 


43-9i±S:IS 


0-70±g:S 


9 c+0.5 
^• J -0.3 


0.97,0.71,0.65 


0.4-0.8 


0.6 


143 


(1.84 ±0.15) 


10" 


■6 


43.911°;!! 




9 o+03 
^• J -0.2 


1.00,0.03,0.08 


0.8-1.6 


1.2 


187 


(1.19 ±0.09) 


10" 


-5 


43.971°;^ 


10 +017 

u - lu -0.19 


9 q+0-2 
0.1 


0.91,0.68,0.88 


1.6-3.2 


2.4 


110 


(1.03 ±0.10) 


10" 


■5 


44.39+°^ 


U - 1O -0.18 


9 O+02 
^• J -0.2 


0.91,0.55,0.44 


3.2-4.8 


4.0 


17 


(4.53 ±1.10) 


■ 10" 


■6 


44 4 q+o:40 


-0 26+ 0,56 


1 q+o.4 
l- y -0.3 


0.95,0.21,0.39 



Parameter values which have been fixed during the fit are labelled by '(*)'. a Units - A44: h% Mpc -3 , L x ,«: 10 44 h^erg s -1 . b 
The three values are probabilities in two 1D-KS test for the distribution, L x , 1D-KS test for the z distribution and the 2D-KS 
test for the (L^,z) space respectively. 



4.3. Evolution of the Space Density 

In this section, we investigate the evolution of the type- 
1 AGN space density in different luminosity classes as a 
function of redshift. The estimator of the space density 
is the N ohs /N mdl . The fit with Eq. El has been made in 
finer redshift shells than in Sect. 14.21 The space densities 
as a function of redshift were calculated in five luminosity 
classes with log L x of 42-43, 43-44, 44-45, 45-46, and >46 
as well as the sum over all luminosities with log L x > 42. 
The resulting curves are shown in Fig. Ufa) for the nom- 
inal calculations. The incompleteness upper bounds have 
also been calculated, but have not been shown here for 
the visibility of the figure. These upper bounds are shown 
in Sect. El (Figs. HU&EJ. Since the Black Hole growth 
function is more closely linked to the emissivity per co- 
moving volume, we also show the emissivity as a function 
of redshift in the same luminosity classes in Fig. EJb). 

Figure EI&) clearly shows a shift of the number den- 
sity peak with luminosity, in the sense that more luminous 
AGNs (QSOs) peak earlier in the history of the universe, 
while the low luminosity ones arise later. Also, there is 
a clear decline of the derived space densities at least for 
luminosities of log L x < 44, even when the optical in- 
completeness upper bounds are taken into account. The 
counting statistics and spectroscopic incompleteness for 
the more luminous AGNs do not allow to determine a 
decline, but do also not exclude it. This issue is further 
discussed in Sections E] and El 

In order to show the behaviour of the luminosity de- 
pendence of the evolution more quantitatively, we have 
also made a maximum-likelihood fit of the evolution curve 
in each of the luminosity bins, with log L x ranges of 42- 
43,43-44,44-45, and 45-46. We used two power law com- 
ponents of the (1 + z) evolution with a cutoff redshift: 



e d (z,L x 



(1 + zjP 1 

e d (z c )[(l + z)/{l + z c )]P 2 (z>z c ) 



p2 f„ ^ , • I 4 ) 



where L xc is the central (logarithmic) luminosity of the 
bin. As was the case for Eq. El the fit depends on the 
shape of the luminosity function (along the luminosity di- 
rection) within the luminosity bin. Again, we have fixed 
the behavior in the luminosity direction using those from 



the LDDE model (Sect. rOjl as a template. We also show 
the normalization: 



A = d$(i x = L xc , z = 0)/dlogi x 



(5) 



The best-fit results are shown in TableEJtogether with the 
K-S probabilities. In the logi x = 45 — 46 bin, the fit for 
P2 was unconstrained. Thus we have fixed the values of P2 
to that from LDDE (see below) for this luminosity bin. 

4.4. Global Representations by Analytical Functions 

It is sometimes useful to provide a simple analytical fit for 
the SXLF over the whole redshift-luminosity range. We 
first used the pure-luminosity evolution (PLE) form, in 
order to enable a comparison with previous work: 



d<P(L x ,z) = d$(£ x /ei(z),0) 
d log L x d log L x 

with the luminosity evolution factor: 



ei(z) 



(1 + zY 1 



(z < z c ) 



e(z c )[(l±z)/(l±z c )f 2 (z>z c ) 



(G) 



(7) 



We again used the smoothed two power law form 
(Eq. El excluding the ^-dependent factor) for the z = 
SXLF. The best-fit PLE parameters are shown in TablcEl 
with the results of the K-S tests. The best-fit PLE model 
is overplotted with the binned SXLF in Figs. El & El as 
dotted lines. It is apparent from the comparison in these 
figures, especially the latter, PLE does not represent the 
behaviour of the low-luminosity (logL x ^ 44), interme- 
diate redshift (0.5 ^ z ^ 1.8) regime, due to the rather 
restrictive nature of the PLE form. 

As a more general analytical form for a refined repre- 
sentation of the SXLF, we have explored the luminosity- 
dependent density evolution form (LDDE) form, origi- 
nally suggested by Schmidt & Green ( 1983) for describing 
optically-selected QSOs: 



d<P(L x ,z) = d$(L x ,0) 
d log L x d log L x 



e<i(z,L x ), 



(8) 



where ed(z,L x ) is the density function normalized to 
z=0. The results from Sect. l4.3l show that the peak number 
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Table 3. The full binned SXLF values 



z 




log 




jyobs 


jymdl 


d# 
d log L x 


b 




(iV f ) 


0.015 - 


0.2 


42.00 


- 42.50 


10 


11.5 


(i-3tr 4 ) ■ 


10" 


5 




0.015 - 


0.2 


42.50 


- 42.75 


15 


15.3 


(i-o+oi) - 


10" 


5 




0.015 - 


0.2 


42.75 


- 43.00 


20 


25.2 


(6.41M) • 


10" 


6 




0.015 - 


0.2 


43.00 


- 43.25 


48 


38.6 


(7.iti:S) • 


10" 


6 




0.015 - 


0.2 


43.25 


- 43.50 


45 


44.9 


(3.41S:?) • 


10" 


6 




0.015 - 


0.2 


43.50 


- 43.75 


38 


45.7 


(i-4t° :5) • 


10" 


6 




0.015 - 


0.2 


43.75 


- 44.00 


36 


34.6 


(6.511:?) • 


10" 


7 




0.015 - 


0.2 


44.00 


- 44.25 


27 


23.7 


(2-4l°'j) • 


10" 


7 




0.015 - 


0.2 


44.25 


- 44.50 


22 


16.9 


(8.5±?:S) • 


10" 


8 




0.015 - 


0.2 


44.50 


- 45.50 


7 


11.2 


(2.i±J:i) • 


10" 


9 





0.2 - 


0.4 


42.25 


- 43.00 


9 


9.6 


0.2 - 


0.4 


43.00 


- 43.25 


12 


11.2 


0.2 - 


0.4 


43.25 


- 43.50 


23 


19.1 


0.2 - 


0.4 


43.50 


- 43.75 


22 


24.1 


0.2 - 


0.4 


43.75 


- 44.00 


25 


25.7 


0.2 - 


0.4 


44.00 


- 44.25 


15 


18.1 


0.2 - 


0.4 


44.25 


- 44.50 


11 


10.6 


0.2 - 


0.4 


44.50 


- 44.75 


11 


7.6 


0.2 - 


0.4 


44.75 


- 45.50 


11 


11.8 


0.2 - 


0.4 


45.50 


- 46.50 





0.4 



2.otg : ?) 
(lo.olij 
7.oll:I 

3.0 ' " 
1.5 
4.7 



0.4 - 


0.8 


42.00 - 


42.75 


18 


17.1 


(i-itS:!: 


■ 10" 


4 


0.4 - 


0.8 


42.75 - 


43.50 


23 


22.0 


(2.3±g : «; 


■ 10" 


5 


0.4 - 


0.8 


43.50 - 


43.75 


13 


15.4 


(5.511:*; 


■ 10" 


6 


0.4 - 


0.8 


43.75 - 


44.00 


21 


23.4 


(2.6l°j; 


■ 10" 


6 


0.4 - 


0.8 


44.00 - 


44.25 


26 


22.7 


(i-3l° :i; 


■ 10" 


6 


0.4 - 


0.8 


44.25 - 


44.50 


13 


16.0 


(2.9±S:g: 


■ 10" 


7 


0.4 - 


0.8 


44.50 - 


44.75 


11 


11.1 


(i-otsi: 


■ 10" 


7 


0.4 - 


0.8 


44.75 - 


45.25 


13 


8.8 


/ 9 9 +0. 7* 
V 4 -^— 0.6/ 


■ 10" 


8 


0.4 - 


0.8 


45.25 - 


46.25 


5 


6.2 


(2.2ll: ) 


■ 10" 


10 


0.4 - 


0.8 


46.25 - 


47.25 





0.2 


< 1.8 ■ 


io- n 





0.8 - 


1.6 


42.00 


- 42.50 


12 


11.0 


(4.211:5: 


■ 10" 


5 


0.8 - 


1.6 


42.50 


- 43.00 


13 


12.4 


(3.611:2; 


• 10" 


5 


0.8 - 


1.6 


43.00 


- 43.50 


14 


17.7 


/ 9 a+0. 8' 
\ ZA -0.6> 


• 10" 


5 


0.8 - 


1.6 


43.50 


- 43.75 


12 


12.9 


In 9+0.8* 
l z ' z -0.6/ 


■ 10" 


5 


0.8 - 


1.6 


43.75 


- 44.00 


20 


16.4 


(2.ol°:i: 


■ io- 


5 


0.8 - 


1.6 


44.00 


- 44.25 


27 


25.5 


(8.211:1; 


■ 10" 


6 


0.8 - 


1.6 


44.25 


- 44.50 


27 


25.0 


In q+0.6n 
( Z -°-0.5, 


■ io- 


6 


0.8 - 


1.6 


44.50 


- 44.75 


26 


25.7 


(7.511-J; 


• io- 


7 


0.8 - 


1.6 


44.75 


- 45.00 


13 


16.7 


(i-5l° :S: 


■ io- 


7 


0.8 - 


1.6 


45.00 


- 45.75 


15 


18.5 


(i-Ho-.l: 


■ io- 


8 


0.8 - 


1.6 


45.75 


- 46.75 


8 


4.4 


(2.2l o: °) 


■ 10" 


10 


0.8 - 


1.6 


46.75 


- 47.75 





0.3 


< 4.0 • 


io- 12 



1.6 - 


3.2 


42.00 


- 42.75 


10 


9.5 


1.6 - 


3.2 


42.75 


- 43.25 


11 


13.1 


1.6 - 


3.2 


43.25 


- 44.00 


23 


22.0 


1.6 - 


3.2 


44.00 


- 44.50 


27 


23.0 


1.6 - 


3.2 


44.50 


- 45.00 


16 


16.9 


1.6 - 


3.2 


45.00 


- 45.25 


12 


8.4 


1.6 - 


3.2 


45.25 


- 46.00 


8 


13.1 


1.6 - 


3.2 


46.00 


- 47.00 





2.9 


1.6 - 


3.2 


47.00 


- 48.00 


3 


0.5 


3.2 - 


4.8 


43.00 


- 44.00 


5 


4.9 


3.2 - 


4.8 


44.00 


- 45.00 


8 


6.9 


3.2 - 


4.8 


45.00 


- 46.00 


1 


3.2 



+0.8\ 
0.6) 

+0.-^ 
0.3) 



-1.2) 

i.8i° :S) 
6.4l?:S) 
2.5i° o: ?) 

< 9.9 • 10 



6 



10" 
10 

io- 

10~ 6 
10~ 6 

io- 7 
io- 7 

10~ 8 

io- 9 

-11 



2.2 4 
1.4 4 
1.41 



+0.9\ 
0.7) 
+0.5\ 
OA) 
+0.3\ 
0.3) 

8.411:?) 

1- 21S1) 

2- 91^) 

9-H 4 3 :i) 

< 1.2 • 10~- 
4.6l|:§) ■ io 



3.2 - 4.8 



46.00 - 47.00 



+2A\ 
1.6) 

i-7l?i) 

1.9 r9.s+?-? s i 



3.8 
3.2 



■ io- 

■ io- 

■ io- 

10" 



(1) 

(6) 



Hasinger et al.: Luminosity-dependent evolution of soft X-ray AGN 



11 



i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i j 

O z = 0.4-0.8 




42 43 44 45 46 47 48 42 43 44 45 46 47 48 42 43 44 45 46 47 48 

Log L x [h-^ergs" 1 ] 

Fig. 4. The soft X-ray luminosity function of the type-1 AGN sample in different redshift shells for the nominal case 
as labeled. The error bars correspond to 68% Poisson errors of the number of AGNs in the bin. The best-fit two power 
law model for the 0.015 < z < 0.2 shell are overplotted in the higher redshift panels for reference. The dotted and 
long-dashed lines give the best-fit PLE and LDDE models discussed in Sect. 14.41 respectively. 



Table 4. Best-fit evolution parameters for each luminosity bin 
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10" 
10" 
10" 
10" 


-a 

-6 
-8 
11 


4.90ti£ 
3-89«;« 

c r 1+0.38 
°- oi -0.37 


0.6512;- 

i 11+0.22 
J " J " L -0.11 
i 7C+0.14 

1 7Q+0-59 


-2.4+^ 

-1.811:2 

-0.4(*) 


0.47,0.77,0.64 
0.55,0.25,0.39 
0.05,0.47,0.07 
0.81,0.98,0.62 



Parameter values which have been fixed during the fit are labelled by '(*)'. a Units - A$: h^ Mpc -3 , L x ,»: 10 44 ft^" 2 erg s _1 . b 
The three values are probabilities in two 1D-KS test for the distribution, L x , 1D-KS test for the z distribution and the 2D-KS 
test for the (L x ,z) space respectively. 
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Fig. 5. (a) The space density of AGNs as a function of redshift in different luminosity classes and the sum over all 
luminosities with log L x > 42. Densities from the PLE and LDDE models (Sect. EOjl are overplotted with solid lines, 
(b) The same as (a), except that the soft X-ray emissivities are plotted instead of number densities. The uppermost 
curve (black) shows the sum of emissivities in all luminosity classes plotted. 



density shifts from z ~ 0.7 at log L x ~ 42.5 to z > 2 
at high luminosities. Based on a similar observation of 
a hard X-ray-selected sample, Ueda et al. 120031 used an 
expression where z c is a simple function of L x : 



e d (z,L x ) 

along with 



(1 + z) pl (z < z c ) 

e d (z c )[(l + z)/(l + z c )}P 2 (z>z c ) 



z c {L x ) 



z c,o(L x / Lx,c) a {L x < L xc ) 
z c,o (L K > L x c ) 



(9) 



(10) 



The results of the analysis in the previous section 
shown in Tabled suggest that considering the dependence 
of pi andp2 on luminosity would still improve the fit. Thus 
we have also included the following for our full LDDE ex- 
pression: 



pl(L x ) = pl M + fa (logL x - 44) 
P 2(L X ) = P 2 44 + (3 2 (logL x -44) 



(11) 
(12) 



The best-fit parameters and the results of the K-S 
tests for the PLE and LDDE models are summarized in 
Table The best-fit PLE and LDDE models are over- 
plotted on Figs. 0] and [5] with dotted and dashed lines 
respectively. A detailed discussion of the comparison of 
model and data is given in Sect. [HJ 



5. An alternate approach using the V max method 

As described in the Introduction, the luminosity function 
derived from survey data binned in luminosity and red- 
shift does not necessarily apply to the centers of the (L x , z) 
bins. This binning bias tends to be especially a problem 
if data are scarce (often at higher redshifts) and gradi- 
ents across bins are large. The previous section describes 
a procedure that corrects the binned space densities to 
first order. 

In this section, we avoid deriving densities from binned 
survey data. Instead, we use the V max values of individual 
RBS sources to derive the zero redshift luminosity func- 
tion. We then derive by iteration an analytical density 
template at various L x values that, together with the zero 
redshift luminosity function, accounts for the observed 
counts and redshifts of the deeper surveys. The end re- 
sult of the procedure is a set of observed values of the lu- 
minosity function that apply to the centers of the (L x , z) 
bins, and that is quite insensitive to the precise template 
employed. A further advantage of employing VJ nax of in- 
dividual sources is that it can be derived for two or more 
selection variables. This allows us to account for the effect 
of a spectroscopic magnitude limit in some of the deeper 
surveys beyond which the redshift is unknown for most of 
the sources. In the first use of Vmax, this feature was used 
to derive the luminosity function of radio quasars from a 
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Table 5. Best-fit parameters for global expressions 





PLE 


LDDE 


z=0 SXLF Parameters 


A44 


1.85 ± .11 10"'' 


2.62 ± .16 10~ Y 


log L x ,* 


43.33 ± .09 


43.94 ± .11 


7i 


0.39 ± .09 


0.87 ± .10 


72 


2.29 ± .09 


2.57 ± .16 


Evolution Parameters 


pl,pl 44 


2.67 ± .15 


4.7 ±.3 


Zc,Z c ,44 


1.70(+.20; -.24) 


1.42 ± .11 


p2,p2 4 4 


-0.2 ± .7 


-1.5 ±.7 


a 




0.21 ±0.04 


l0gi x ,c 




44.67* 


ft 




0.7 ±0.3 






0.6 ±0.8 


K-S Probabilities (L x ,z,2D) 


0.015 < z < 4.8, lo 


gL x > 42 0.10,0.10,0.05 


0.86,0.65,0.36 



Parameter values which have been fixed during the fit are labeled by '(*)'. a Units - A: h% Mpc -3 , L x ,„: 10 44 /i^erg s" 1 . b 
The three values are probabilities in two 1D-KS test for the distribution, L x , 1D-KS test for the z distribution and the 2D-KS 
test for the (£ x ,z) space respectively. 



sample in which only the optically brightest objects had 
redshifts (Schmidt 



5.1. Using V max to derive the luminosity function 

The derivation of a luminosity function from objects in a 
well defined sample usually involves binning the observa- 
tions in redshift and luminosity. If we make the bins in 
luminosity so small that each contains only one or zero 
objects then the luminosity function is composed of con- 
tributions from each of the individual sample objects. In 
the limit, each of the objects contributes to the luminosity 
function a delta function of amplitude 1/Vmax at the ob- 
ject's luminosity L, where Kn ax is the co-moving, density- 
weighted volume over which the object can be observed 
within the sample limits in flux and solid angle. This lu- 
minosity function will reproduce the source counts of the 
input sample exactly. 

We write the luminosity function as 



$(L x ,z) = $(L x ,0)p(z,L x ) 



(13) 



where p(z, L x ) is the space density, or density evolu- 
tion, normalized to z=0. We approximate p(z, L x ) by an 
analytical density template. For the redshift dependence 
we use at low redshifts a power law of (1 ± z); at higher 
redshifts we adopt the shape of the density function used 
by Schmidt et al. (1995) for optically selected quasars: 



Ptend» = (1 + *) m 0<Z<Z C (14) 

Ptom(z) = (1 + Z c ) m Z c < Z < Z d (15) 

PtcmW = (l±2 c )" l 10' £ ( Z -^ Z>Z d (16) 

We use the RBS which contains 205 AGN-1 in the 
log L x range 42-46, mostly at low redshifts, to derive the 



S o 




lim.x 

Fig. 6. The X-ray redshift limit versus the optical redshift 
limit for 205 RBS AGN-1 sources. In this illustration the 
X-ray flux limit is 10~ 16 cgs and the magnitude limit R = 
24, similar to the case of the CDF-N. 



zero redshift luminosity function $(L X ,0). The main ex- 
cercise then is to derive the values of the template pa- 
rameters (m, z c , Zd, k) by fitting to the flux and redshift 
distributions in the deeper X-ray surveys. This will allow 
the direct derivation of the luminosity function at the cen- 
ter of bins of luminosity and redshift, as outlined in the 
steps below. 

1. Start the iteration by assuming initial values of the 
template parameters (m, z c , z d , k) as a function of L x ; 

2. The zero redshift luminosity function is the sum of 
delta functions for the assembly of RBS sources 



*(ix,o) = E-=i(i« s ,)^i-^ 



(17) 
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Fig. 7. The z=0 luminosity function $(L X , 0) derived from 
205 RBS AGN-1 sources using the V max method (red filled 
symbols), compared with the binned estimate in the red- 
shift shell z=0.015-0.2, corrected to z=0. 



where j is the accessible density-weighted volume 
of RBS source i in the RBS, based on its solid angle 
and flux limit distribution; 

Next, predict expected numbers and redshifts in the 
deeper surveys. For the part of the luminosity function 
based on RBS source i, the expected number in survey 
sur based on the density template is 



(18) 



where V sur . 

v ^ ^ 'max.? 



5. 



The values of V™Z , are based on the maximum accessi- 

III cl X , L 

ble redshifts i of RBS source i in survey sur. For the 
SA-N, NEPS, RIXOS, RMS and CDF-S surveys, where 
missing redshifts are not correlated with X-ray or optical 
flux, we assume that z^[ x i equals the X-ray redshift limit 
ziim.x- In accounting for the missing redshifts, we assume 
that their distribution is similar to that of the observed 
redshifts and use an effective solid angle appropriately re- 
duced from the geometric solid angle. 

For the RDS/XMM and CDF-N surveys, where red- 
shifts are missing primarily for optically very faint objects, 
we proceed as follows. In this case, we employ only objects 
brighter than a magnitude limit mn m (in practice R < 24) 
and use the full geometric survey solid angle. In predicting 
the number of expected objects in survey sur, this intro- 
duces a second limiting redshift, zii mj0 p tj i that will depend 
on the optical luminosity of RBS source i contributing to 
the luminosity function. The relevant limit z^[ x i for the 
derivation of the expected counts is the smaller of zn m . x ,i 
and ziim,opt,i- 

The determination of zii m! opt,i is based on the mag- 
nitudes of the RBS objects. Among the 205 RBS AGN- 
1, Salvato 1)2002(1 has carried out optical photometry of 
a redshift-limited sample of 89 sources, deriving magni- 
tudes for the nucleus, disk and spheroid. Comparison of 
the total magnitudes given by Salvato with the RBS mag- 
nitudes shows < -Rrbs — Rs&\ >= —0.5 with a dispersion 
of 0.9 mag. We apply this correction to the 116 remain- 
ing RBS sources. Based on the systematics of the Salvato 
magnitudes, we find that on the average, the absolute R- 
magnitude of the galaxy light is -21.3, that the disk light 
accounts for 56% of the galaxy light, and that for low- 
luminosity objects the nucleus contributes no less that 
25% of the total light. This allows the, admittedly un- 
certain, derivation of the nuclear, disk and spheroid mag- 
nitudes of the remaining 116 RBS objects. We use the 
following values for the spectral index a = d\ogS / d\ogv 
of the nuclear, disk and spheroid components, estimated 
from spectral energy distributions given by Kinney (Q-996 ): 
c^nuc = —0.5, adis = —1-0, a sp h = —3.0. These values al- 
low evaluation of the K-correction required to derive the 
redshift zn m ,opM a ^ which each of the 205 RBS sources 
reaches the spectroscopic magnitude limit i?ii m = 24. We 
illustrate in Figure the effect of the optical limit on the 
redshift limits for the CDF-N: for objects below the di- 
agonal line, the redshift limit is that set by the optical 
limit. 



(z) is the accessible density-weighted vol- 
ume of RBS source i over the redshift range z, z + Az 
based on the solid angle and flux limit distributions in 
survey sur; 

In order to compare with the observed numbers, 
we use four luminosity classes with AlogL x = 1.0 
centered on logL x = 42.5,43.5,44.5,45.5 and red- 
shift shells of Alogz = 0.1 centered on logz = 
—0.95, —0.85, 0.75. Total binned numbers predicted 
in all the non-RBS surveys are iVt om (L x , z); 
Next, use the observed number iV bs(£ x , z) in each 
(L x ,z) bin to derive the space density or luminosity 
function at the center of each luminosity bin by scal- 
ing the template value of the luminosity function: 

5.3. The zero redshift X-ray luminosity function 

$ obs (i x ,z) = $(L x ,0) Ptcm (z,L K )N ohs (L x ,z)/N tcm (L x ,z)£\Q) zeTO redshift i um i nos i ty function $(L X ,0) is derived 



6. Comparison of $ ohs (z,L x ) and ^(L x ,0)p tcm (z, L x ) 
serves as a guide for the next iteration of the template 
parameters (m, z c , Zd, fc), starting at step 2 above. 



from the RBS sources as outlined in Sect 15.11 Eq. El 
Since the accessible volume V|]^ S j is density- weighted, this 
derivation requires information about the density tem- 
plate, discussed in the next section. At the low redshifts of 
the RBS objects (typically 0.1), this involves only the m 
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Table 6. Parameters characterizing the density template. 



parameter of the low redshift (l±z) m density variation. A 
binned version of <1>(L X , 0), created by summing the delta 
functions in bins of AlogL x = 0.2, is shown in Fig.[7| 

5.4. The analytical density template p(z,L x ) 

The density template defined in Eas. I15ITTH1 with param- 
eters (m, z c , Zd,k) is an analytical approximation of the 
density function p(z,L x ). The derivation of the template 
parameters was carried out as follows. As described in 
Sect. 15.11 it is a procedure of trial and error. For each of 
the four luminosity classes, in iterating the value of m, we 
minimize < N b s (L x , z) — N tcm {L x , z) > for 0.5 < z < z c , 
since at low redshifts the luminosity function is firmly an- 
chored by the RBS objects. The values of z c , z^ and k are 
derived by minimizing < iV bs(£x, z) — Ar tem (L x , z) > for 
z > z c , by adjusting k for the lower luminosity classes to 
fit the total observed AGN-1 in the CDF-S and CDF- 
N (regardless of availability of redshifts or identification), 
and for the higher luminosity classes by adopting parame- 
ter values for optically selected quasars (see next section) . 
This procedure makes it difficult to evaluate the statistical 
significance of the templates. We will show below that the 
accuracy of the templates actually has a negligible effect 
on the values of <E> bs(£x, z). 

The resulting values of (m, z c , z<i, k) for the four lumi- 
nosity classes are given in Table [SJ The errors on for the 
low-redshift power index m = rflogp/dlog(l ± z) have 
been estimated using the number of observed objects in 
the redshift shells 0.5 < z < z c in comparison with the 
RBS objects. From this analysis and the results shown in 
Table 01 it appears that m exhibits a significant increase 
with L x for logL x > 43.0. For these luminosities, we em- 
ploy a quadratic interpolation of m with logL x through 
the values given in the table. For logL x = 42 — 43 we 
assume m to be constant. 

5.5. The X-ray luminosity function & bs{L x , z) 

The observed values of the luminosity function Q hs{Lx, z) 
are obtained by scaling the template luminosity function 
by a factor N \, S (L X , z)/N tem {L Xl z), see Eq.^D These val- 
ues of $obs(-kx, z) are plotted versus the redshift for each 
of the four luminosity classes in Figs.EHSI The ±lcr error 
bars are based on the numbers N tem (L x , z) predicted by 
the template. The figures also show the template lumi- 
nosity function resulting from the product of the zero red- 
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log L x = 44-45 




z 

Fig. 10. The luminosity function versus redshift z, derived 
using a template having constant space density at large 
redshift. The observed numbers at z = 4.47, 5.62, 7.08 are 
1, 0, 0, respectively. This plot gives clear evidence that 
high-luminosity X-ray AGNs type 1 show a decline in 
space density beyond z ~ 4. 

shift luminosity function $(L X ,0) and the template den- 
sity function p tcm (z,L x ). 

The ±lcr error bars do not include any contribu- 
tion reflecting the error of the template density function 
Ptem(z, L x ). We explore the effect of the template by red- 
eriving the predicted number Nt em (L x , z) on the extreme 
assumption that m = and k = 0, i.e that there is 
zero density evolution. We find that the observed densities 
log $ bs(£x, z) increase by only 0.00 — 0.06. Since the tem- 
plate errors must be much smaller than assumed in this 
extreme example, their effect on the error of <& b s (Lx, z) 
is negligible. The main result of this section is the set of 
®obs(Lx, z) values plotted in Figs.\^Q the templates serve 
primarily to eliminate the uncertainty related to binning. 

For log L x = 42 — 43 and 43 — 44, the density rises by 
an order of magnitude to z ~ 0.7 and 1.2, respectively, 
and then declines steadily (Fig. [SJ . As shown in the next 
section, the peak at z ~ 0.7 for logi x = 42 — 43 is little 
affected by redshift spikes caused by large-scale structure. 

The density distribution for AGN-1 with log£ x =44- 
45 is documented to z ~ 4 (Fig. EJ). Since in this case, the 
evidence for a decline in density at high redshift was not 
initially clear, we adopted for the parameters (z c ,Zd,k) 
of the template above z = 1.7 those found for high- 
luminosity optical quasars by Schmidt et al. (1995). It 
appears from Fig. that the X-ray data are consistent 
with the adopted shape. In order to further investigate 
whether the space density declines significantly at high 
redshift, we explore a test template in which the density 
does not decline at all, i.e. remains flat at high redshifts 
(k = 0). The results are shown in Fig.^3 The error bars 
are much reduced from those in Fig. reflecting the fact 
that our errors are based on the predicted numbers, which 
are larger for a flat evolution function at high redshift. 



logL x 42.5 43.5 44.5 45.5 ft 

C3 

m 4.0 ±0.7 3.4 ±0.5 5.0 ± 0.2 7.1 ± 1.0 ~ 

2 C 0.7 1.2 1.7 1.7 J 

z d 0.7 1.2 2.7 2.7 o 

k -0.32 -0.32 -0.43 -0.43 Z 
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Fig. 8. The luminosity function as a function of rcdshift z. The z = values are based on 21 RBS sources for 
log L x = 42 — 43, 111 sources for log L x = 43 — 44. The points are based on numbers observed in the non-RBS surveys. 
The thin line represents the density template (see text for details). 




Fig. 9. The luminosity function as a function of redshift z. The z — value are based on 67 RBS sources for 
logi x = 44 — 45, 6 sources for logL x = 45 — 46. The points are based on numbers observed in the non-RBS surveys. 
The thin line represents the density template; the parameters of the template for z > 1.7 were adopted from those 
found for high-luminosity optical quasars (see text for details). 



The bins at z = 4.47,5.62,7.07 have observed/predicted 
numbers of 1/7.9, 0/7.2, and 0/19.3, respectively. Limiting 
the case to z < 6.3, for which the observation of redshifts 
in high-luminosity quasars should be no problem, we have 
for z > 4 one observed object versus 15 expected. The 
Poisson probability for such an occurence barring system- 
atic effects is 3 x 10~ 7 , constituting strong evidence that 
for logL x = 44 — 45 the space density declines beyond 
z ~ 4. 

The low-redshift density parameter m for log L x = 45— 
46 continues the trend of an increasing m for larger L x . At 
large redshift, the sparse data do not give any information 
about the density beyond z ~ 3 (see Fig. EJ). 

6. Discussion 

In the present paper we have used two different meth- 
ods to derive the AGN-1 X-ray luminosity function and 
its evolution. Detailed descriptions of the two methods 



are given in Sects. FTTI and PTTI Conceptually, the binned 
method derives a first order luminosity function by divid- 
ing the numbers N a i, s (L x , z) observed in the input surveys 
by the appropriate volumes. An analytical representation 
of the luminosity function is used to predict the num- 
bers N m di(L x , z) expected in these surveys. The luminos- 
ity function is then corrected by the factor N b s /N m di- 

In the V max method, the RBS is used to derive the 
zero redshift luminosity function. An analytical density 
template is used to predict the numbers N tcm (L x , z) ex- 
pected in the deeper surveys. In this process, the effect of 
a spectroscopic magnitude limit on the V max values is in- 
cluded for deep surveys where this limit applies. Once the 
template predictions are close to the observed numbers 
N Q hs(L x , z) in the deeper surveys, the luminosity function 
at the center of each (L x , z) bin is derived by multiplying 
the template luminosity function by N b s /N te m- 

Using the optical magnitudes and spectra of the RBS 
sources in deriving the redshifts at which they would reach 



Hasinger et al.: Luminosity-dependent evolution of soft X-ray AGN 



17 



the spectroscopic limit R = 24 introduces uncertainties. 
The magnitudes of the 116 RBS sources not studied by 
Salvato (2002J) are quite poor and so is their separation 
in nucleus, disk and spheroid components. The assumed 
spectra of these components are schematic. Obtaining 
spectral energy distributions for all the RBS sources down 
to the far UV would allow deriving their flux directly at 
or near the rest wavelengths corresponding to the R mag- 
nitude at the limiting redshifts. 

It is reassuring that the general properties and abso- 
lute values of the space density are very similar in the two 
different derivations in Sections 01 and FigurelTTlshows a 
direct comparison between the binned and VJ nax determi- 
nations of the space density, which agree very well within 
statistical errors. 

We use overall fits to the luminosity function and its 
cosmological evolution in order to compare to previous 
work and to enable theoretical calculations with simpli- 
fied analytical forms. For this purpose, we use two func- 
tional forms, i.e. a pure luminosity evolution (PLE) and 
a luminosity-dependent density evolution (LDDE) model. 
Assuming PLE, the luminosity evolution index pi w 2.7 
and the cutoff rcdshift z c « 1.7 are in rough agreement 
with previous results from type-1 AGNs in ROSAT or 
Einstein/ROSAT combined surveys (e.g. paper I; Jones 
et a UTMTI Page et al lTMTl Boyle et al. lT5M|> . 

The results of the K-S test for the PLE fit over 
the whole sample are marginally acceptable, with ~ 5% 
chance of obtaining a 2D K-S value larger than observed. 
However, this is caused by the sheer number of AGNs in 
the part of the z-L x space where PLE is still a good de- 
scription, which dominates the overall statistics. Figure 
clearly shows that PLE fails to reproduce the behavior 
of the SXLF around log£ x ^ 44 at 0.4 <, z <, 1.7. The 
luminosity bin 42 < logL x < 43.5 in the 0.4 < z < 0.8 
shell contains 41 objects where PLE predicts 19 and the 
adjacent regime of 43.0 < logL x < 44.0 at 0.8 < z < 1.2 
contains 46 objects where PLE predicts 25, so that the 
overall excess corresponds to more than 4.5cx. In construct- 
ing the LDDE form, we tried to fully represent the SXLF 
from our data with unprecedented redshift and luminos- 
ity coverage. The overall 2D K-S acceptance of LDDE 
has been improved to 36%. The only location where the 
LDDE model still deviates from the data significantly is 
the very end of our sample, the 42 < log < 43 bin at 
0.015 < z < 0.2, where LDDE predicts 66 objects while 
we observed 45 objects (a ~ 3a deviation). 

Even though our sample is a soft X-ray-selected type-1 
AGN sample, the overall behaviour of our XLF is similar 
to that obtained by Ueda et al. (2003) for the intrinsic 
(de-absorbed) luminosity function of hard X-ray selected 
obscured and unobscured AGN. To make this compari- 
son, we have refitted our sample with the LDDE model 
where j3\ and /?2 are fixed to zero. This is exactly the 
same as the function form which Ueda et al. I|20031 used 
to describe their intrinsic HXLF. All the z = XLF pa- 
rameters and evolution parameters are remarkably close 
between our SXLF and HXLF, except the global normal- 



ization. The HXLF normalization is found to be about five 
times larger than that of our SXLF, after adjustments for 
the differences in energy bands and the difference in the 
luminosities at which normalizations are evaluated. This 
factor probably accounts for the absorbed objects missing 
in the SXLF. However, the Ueda et al. sample, contain- 
ing about 250 AGN, is limited to lower luminosities and 
lower redshifts than our sample of 1000 objects, so that its 
statistical quality and limited sensitivity range were not 
sufficient to constrain the decline of the space density at 
high redshifts, which has been measured significantly in 
our sample for the first time. 

Very recently, Barger et al. , 2005jhave presented X-ray 
luminosity function analyses both in the hard and soft X- 
ray bands, based on the CDF-N, CDF-S, CLASXS and 
ASCA surveys. Again, their results are in good agreement 
with the soft XLF discussed here and the hard XLF pre- 
sented by Ueda et al., however, they still suffer from sub- 
stantial identification incompleteness. Also, their results 
on broad-line AGN are not directly comparable to our 
type-1 AGN sample, because they only include the op- 
tically classified type-1 AGN and thus miss most of the 
low-luminosity unabsorbed AGN-1 we are concerned with 
in this paper. A critical comparison of our XLF with those 
of Ueda et al.. l2003l and Barger et al.. l2005l will be the topic 
of a future paper. 

The evidence for a peak in the evolution function is 
quite strong at z ~ 0.7 for L x = 42 — 43 (see Figs. CJa) 
and El 

The faintest end of the sample depends on the small 
field of Chandra Deep Field-South and there is some con- 
cern on the effects of the large-scale structure and cos- 
mic variance associated with it. In particular, Gilli et al. 
2003} found two redshift spikes, one at z ~ 0.67 and the 
other at z ~ 0.73. One may wonder whether the z=0.7 
peak is caused by these redshift spikes (see also Gilli et 
al., 2003). These do, however, not affect our SXLF esti- 
mates significantly. In the analysis in Sect. 0] only six out 
of 41 AGNs (15%) in our sample in the 0.4 < z < 0.8 
bin in the range 42 < logL x < 43.5 are in these spikes 
(0.664 < z < 0.685 and 0.725 < z < 0.742) and this is the 
only regime where the spikes give a non-negligible contri- 
bution. In the analysis in Sect. [5] the z=0.71 bin contains 
5 objects from the two spikes, for a total of 5 out of 15 ob- 
served objects. Disregarding these objects would decrease 
the derived density by dex 0.22 in Fig. actually lead- 
ing to better agreement with the template densities. We 
therefore conclude that cosmic variance is not significantly 
affecting our results on the evolution of the space density. 

We show in Figure ITSl the AGN-1 space density as a 
function of cosmic time. We see dramatic changes with L x . 
For declining L x as we move from high-luminosity AGN 
or quasars to Seyfert galaxies, the main formation of the 
objects occurred at later cosmic times. For logL x > 44 
the density curve is similar to that for quasars. It is an 
intriguing question whether the observed dependence of 
m on L x is accompanied by a corresponding dependence 
on L opt . At logL x < 44 the AGN-1 are mostly Seyfert 
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galaxies, for which there is no comparable optical evidence 
about their density curve. 



These new results paint a dramatically different evo- 
lutionary picture for low-luminosity AGN compared to 
the high-luminosity QSOs. Obviously, the rare, high- 
luminosity objects can form and feed very efficiently rather 
early in the universe. Their space density declines by more 
than two orders of magnitude at redshifts below z=2. The 
bulk of the AGN, however, has to wait much longer to 
grow with and shows a decline of space density by less 
than a factor of 10 at redshifts below one. The late evo- 
lution of the low-luminosity Seyfert population is very 
similar to that which is required to fit the Mid-infrared 
source counts and background (Franceschini et al., 2002) 



and also the bulk of the star formation in the Universe 
(Madau et al.. I1998j) . while the rapid evolution of pow- 
erful QSOs traces more the merging history of spheroid 
formation (Franceschini et al. 1999). 



This kind of anti-hierarchical Black Hole growth sce- 
nario is not predicted in most of the semi-analytic models 
based on Cold Dark Matter structure formation models 
(e.g. Kauffmann & Haehnelt MEM Wyithe & LoebHHSI). 
This could indicate two modes of accretion and black hole 
growth with radically different accretion efficiency (see 
e.g. Duschl & Strittmatter |2002J) . A self-consistent model 
of the black hole growth which can simultaneously ex- 
plain the anti-hierarchical X-ray space density evolution 
and the local black hole mass function derived from the 
•^bh - cr relation assuming two radically different modes of 
accretion has recently been presented by Merloni l|2004[) . 

Finally, we compare the space density of soft X-ray 
selected QSOs from our sample to the one of optically- 
selected QSOs at the most luminous end. The comparison 
is plotted in Fig. 1131 The z < 2 number density curve for 
optically selected QSOs (M^ < —26.0) is from the combi- 
nation of the 2dF and 6dF QSO redshift surveys by Croom 
et al. (2004). The z > 2.7 number densities from Schmidt, 
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Fig. 13. Comparison of the space density of luminous QSOs between optically selected and soft X-ray selected samples. 
The X-ray number densities are plotted for the luminosity class log L x — 44 — 45, both for the binned and T^nax analysis 
with the same symbols as in Fig.^J The dashed lines in represent the one sigma range for M&. < —26.0 from Croom 
et al 120041 multiplied by a factor of 40 to match the X-ray space density at z=1.5. The triangles at z > 2.7 with la 
errors are from Schmidt, Schneider, & Gunn (1995) (SSG95) and Fan et al. (2001) after a cosmology conversion (see 
text) and a scaling by a factor of 40 to match with the soft X-ray density at z ~ 2.7. As discussed in this paper, both 
the rise and the decline of the space density change with L x and therefore this comparison can only be illustrative. 



Schneider & Gunn (|1995[) and Fan et al. H2001[) have been 
originally given for hi = 5/7, Sl m = 1,Oa = 0. Their 
data points have been converted to our default cosmology 
and the Mb threshold has been re-calculated with an as- 
sumed spectral index of a = —0.79 (/„ oc ^ Q °), following 
e.g. Vignali et al. (2003- The plotted curve from Schmidt, 
Schneider & Gunn / Fan et al. is for Mb < —26.47 un- 
der these new assumptions. A small correction of den- 
sities due to the cosmology conversion causing redshift- 
dependent luminosity thresholds has also been made, as- 
suming f/'l'. f/log/.;. x /.,. ; '' (Fan et al. 200 Q. The space 
density for the soft X-ray QSOs for the luminosity class 
44 < log L x < 45 has been plotted both for the binned 
and Vmax determination. The Croom et al., (2004} space 
density had to be scaled up by a factor of 16 in order to 
match the X-ray density at z ~ 2. The Schmidt, Schneider 
& Gunn / Fan et al. data points have been scaled by a 
factor of 40 to match the soft X-ray data at z — 2.7 in 



the plot. There is relatively little difference in the density 
functions between the X-ray and optical QSO samples, 
although we have to keep in mind, that both the rise and 
the decline of the space density is varying with X-ray lu- 
minosity, so that this comparison can only be illustrative 
until larger samples of high-redshift X-ray selected QSOs 
will be available. 

Very recently, Wall et al. 120051 have presented an up- 
date of the space density evolution of the Parkes quarter- 
Jansky sample of flat-spectrum radio sources. They basi- 
cally confirm the rise and fall of the QSO population as 
now seen both in the optical and X-ray QSO populations. 

7. Conclusions and outlook 

We have merged the Chandra and XMM-Newton deep 
survey data with the whole body of previously identified 
RO SAT AGN samples. We have selected only the type-1 
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AGN in all samples and treated only the detections and 
X-ray fluxes in the 0.5-2 keV band. The different sam- 
ples cover an unprecedented five orders of magnitude in 
flux limit and six orders of magnitude in survey solid an- 
gle between the ROSAT Bright and serendipitous surveys, 
the XMM-Newton Lockman Hole survey and the Chandra 
Deep Surveys. The sample comprises 944 identified AGN- 
1 objects and only 57 unidentified sources, which could be 
AGN-1, i.e. roundabout 1000 objects. The luminosity- 
redshift diagram is almost homegeneously filled with our 
sample objects. With this sample we arrive at the follow- 
ing conclusions: 

1. The new Chandra and XMM-Newton sources are pre- 
dominantly Seyfert galaxies at a median luminosity of 
~ 10 43 erg s _1 and a median redshift around 0.7 and 
push the determination of the X-ray luminosity and 
space density functions into so far unexplored param- 
eter ranges of redshift and luminosity. 

2. AGN-1 are by far the largest contributors to the soft 
X-ray selected samples. Their evolutional properties 
are responsible for the break in the total X-ray sources 
counts in the 0.5-2 keV band. 

3. The soft X-ray luminosity function of AGN shows a 
clear change of shape as a function of redshift, con- 
firming earlier reports of luminosity-dependent density 
evolution for optical quasars and X-ray AGN. 

4. The space density function changes significantly for 
different luminosity classes. It shows a strong positive 
evolution, i.e. a density increase at low redshifts up to a 
certain redshift and then a flattening. The redshift, at 
which the evolution peaks, changes considerably with 
X-ray luminosity, from z « 0.5 — 0.7 for luminosities 
logi x = 42-43 erg s" 1 to z « 2 for logi x = 45-46 
erg s _1 . 

5. The amount of density evolution from redshift zero to 
the maximum space density is also a strong function 
of X-ray luminosity. The change is more than a factor 



of 100 at high luminosities, similar to what has been 
observed for optically selected QSOs, but it is less than 
a factor of 10 for low X-ray luminosities. 
6. For the first time, we find a clear decline of the space 
density of X-ray selected AGN towards high redshift, 
using a rigorous treatment of optical incompleteness 
and the corresponding survey volume. The decline is 
observed clearly for X-ray luminosities in the range 
\ogL rmx = 42—45 erg s _1 , while at higher luminosities 
the survey volume at high-redshift is still too small to 
obtain meaningful densities. 

In the future, X-ray surveys which are both wide and 
deep are necessary, in order to provide enough volume for 
a better measurement of the space density function of the 
rare high-luminosity AGN at large redshifts. Several new 
surveys towards this goal are already underway, e.g. the 
Chandra Multiwavelength Project (Champ) (Silverman et 
al.. 12003)) . the XMM-Newton AXIS project (Barcons et al., 
2002), the Chandra Large Area Synoptic X-Ray Survey 
(CLASXS) (Yang et al., 127)041 Steffen et al., rSOOll) . the 
Extended Chandra Deep Field South (PI: W.N. Brandt) 
or the XMM-Newton COSMOS Field (PI: G. Hasinger), 
which together should enrich the sample of z > 4 objects 
by about an order of magnitude. Ultimately, new X-ray 
Dark Energy missions, aiming to survey large solid angles 
on the sky to considerable depth could provide a factor of 
100-1000 increase in AGN sample size. 
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